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SUMMARY 

2 . 93*5 

A Fortran IV computer program for calculating the flow fields in three- 
dimensional axisymmetric or two-dimensional inlets has been written. The pro- 
gram has been written to handle inlets designed with smooth compression 
surfaces and for which the attached bow shock falls outside the cowl lip. The 
method of characteristics has been used to calculate a uniform field of points 
and at each of these points the total pressure, Mach number, local flow angle, 
and static pressure ratio are printed. The numerical procedures used are fully 
described and a test case is presented. 


INTRODUCTION 



The mathematical tools for calculating flow fields in supersonic inlets 
have been available for many years. The complexity of the mathematical proce- 
dures, however, has been a major obstacle in effectively and rapidly designing 
inlets other than those with simple two-dimensional compression surfaces. A 
Fortran IV computer program employing the method of characteristics for a per- 
fect gas has therefore been written to assist in the design of three-dimensional 
axisymmetric or two-dimensional inlets. The program is limited in application 
to designs in which the bow shock wave does not intersect the cowl so that 
internal shock wave intersections do not occur (see fig. 1). In addition, no 
viscous effects are considered. Within these limitations the flow fields of 
inlets with theoretical efficiencies up to 100 percent (lOO-percent total- 
pressure recovery) can be described. 

This program has been made available to a number of organizations and the 
purpose of this report is to aid those organizations and any others using this 
program. The basic equations are presented and the program is described fully 
in this report. This description includes the program listing, program usage, 
flow charts, and a sample case. 

The input to the program consists of the surface contours, the free -stream 
Mach number, and other pertinent parameters which are described in appendix A. 
The output consists of a uniform field of points, at each of which the total 
pressure, Mach number, local flow angle, and static pressure ratio are printed. 

This program will be distributed in card form by the Ames Research Center 
upon written request. This distribution also includes complete sample cases, 
program listings, and a complete set of output. 
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SYMBOLS 


a 


M 

P 

r 

S 

V 

W 


x 


7 

b 

0 


1 

2 
l 
t 

U 

o o 


speed of sound 

dimensionless incremental distance along a character- 
istic 

Mach number 
pressure 
radial distance 

entropy divided by universal gas constant 
velocity 

ratio of local velocity to stagnation speed of sound, 
V/a-fc 

axial distance 
ratio of specific heats 
stream angle, radians 
Mach angle, radians 
shock-wave angle, radians 


Program label 

4 

C$>R1, C|)R2 

EM ; EMIN 
PI, p 3; P 1 *- 
R 
S 

W 

X 

GAMMA 

DEL, DELTA 
U 

THETA 


Subscripts 

calculation along a first -family characteristic line 

calculation along a second -family characteristic line 

local conditions 

stagnation conditions 

upstream conditions for a shock wave 

free -stream conditions 
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PROGRAM DESCRIPTION 


The method of characteristics is a standard procedure used in the study 
of supersonic flow fields. It is fully described in references 1 and 2. The 
basic equations used here are essentially those developed in reference 3 but 
with modifications made by Mr. Leroy Presley of Ames Research Center. The 
numerical techniques used in this program are described below. Appendix A 
contains the program usage , a description of the required input, a list of 
error messages, the program listing, and a test case. Because flow diagrams 
are often helpful in remedying program difficulties, a complete set of diagrams 
is presented in appendix B. 


Basic Equations 

The basic equation used in the program is the compatibility equation 


where 


dS = ±A dW + C ± D dS 


(la) 


W = 


V_ 

a t 


A 





n c . 

C = — sin p sin S 


- 1 . n 
D = — sm 2 g 


The upper signs in equation (la) are used along the first-family characteristic 
lines defined by equation (lb). The lower signs in equation (la) are used 
along the second -family characteristic lines defined by equation (lc). 


= tan (p + &) (lb) 

= -tan(p - S) ( lc ) 

The equation for the two-dimensional characteristics program is obtained by 
setting the term C equal to zero. J 
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Basic Computational Techniques 


In order to facilitate the computation, the flow field behind the bow. 
shock wave is broken into several regions bounded by shock waves as shown in 
figure 1. A second -family region exists behind a down-shock (regions 2 and 4 
in fig. 1) and a first -family region behind an up -shock (regions 1 and. 3 in 
fig. 1). If the signs in equation (l) are reversed, the same computation 
schemes may be used for both the first- and second -family regions. Thus, in 
the equations to follow, a double sign implies that the upper. sign is to be 
used in first-family regions and the lower sign in second-family regions. 

In each region successive rays are computed from a surface to a shock, 
wave until the shock wave intersects a surface or falls outside the cowl lip, 
as in the case of the first region. As soon as the intersection occurs a new 
region is started and the previous region continued only in the. area in which 
it is needed, thereby eliminating unnecessary calculations. This is described 
more fully in the section on Shock Point Calculation. 

A numbering scheme has been set up such that each point in the flow field 
is defined uniquely by a region number, i, a ray number, j, and a point num- 
ber, k. (See fig. 2.) In this figure, j ' = j - 1, where j is the. current 
ray number, and j ' is the previous ray number. The subscript k is the 
current point number in the current ray. The subscript k' = k - 1 if j is 
odd and k' = k if j is even. In a first -family region, k' is a first - 
family point; that is, it is located on the first-family characteristic to the 
point k, and k' + 1 is located on the second -family characteristic. In a 
second-family region the k' and k' + 1 points are interchanged. Each odd- 
numbered ray contains a body point. Each ray contains a shock point until 
such time as a body -shock or cowl-shock intersection occurs in the region. 
Three adjacent vertical computation rays are shown in figure 1 with their con- 
necting characteristic lines . 


Field Point Calculation 

In computing a flow-field point, the geometrical location of the point is 
at the intersection of the first- and second -family characteristic lines. The 
x coordinate is found from equation (2) : 

r k' + i ~ V 1 [ x k' tan ^ 1 S) k' + X k'+i tan( ^ + 8) k' + J (2) 

k ± [tan(p ± 5) k , + tan(p + & ) k i +1 ] 

The r coordinate is found by equation (3): 

r k = r k . ± ( x k ' x k’) tan ^ 1 & )k' ^ 
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The distance along the characteristic lines is found from equations (h): 



(£) .(£) . 

r k “ r k' 



r /i 

r k , sin(p ± B) k , j 



} (4) 


( z) = (£) 

r k " r k'+i 

V r /2 V r /k'+i-^k 

r k«+i sin(p + 6) k , +i 


In order to find the stream angle, entropy, and velocity, the compatibility 
equation, equation (1), is put into finite difference form, equations (5): 


6 k - 5 k , = ±A k ,(W k - W k >) + C k . ± D k .(S k - S k ,) 

5 k " & k'+i = + A k’+i( w k " w k '+i) 4 c k'+i + D k'+i( s k 


> 


®k* +i) 


( 5 ) 


This set of two equations in three unknowns may be solved for £>> in terms 
of S k , equation (6): K 



( 6 ) 


By assuming a linear variation in entropy along a normal to the streamline 
through the point k, we may obtain an additional equation. Equations (7) are 
derived geometrically and are shown schematically in figure 3 . 


S k " S k' + a + b ( S k'+i " s k'^) 

a = (r) sin ^ 1 s )k' + & k] > 

b = ^ sin[(p + 6) k , +1 ± 5 k ] j 


( 7 ) 


5 


An initial value of & may be computed by assuming no entropy loss along the 
characteristic lines; that is, the terms (S k - S k « ) and (S k - S k « n) inequa- 
tion ( 6 ) are zero. Equations ( 6 ) and ( 7 ) are then solved iteratively until 
successive values of & converge. The velocity W k may be obtained from 
equations (5) . The remaining properties are then computed . from standard rela- 
tionships found in reference 4 and are presented in equations ( 8 ) through (.13) 


Mv = 



l)W k 2 


(8) 


= sin_1 (j£ 

^*t k 

- — ~ e 




P co V t k / V t*/ 




(9) 



(10) 


7/(7- 1) 


l)W k 2 


(11) 


7/(7-!) 


• DMoo 2 


(12) 

(>) 


(13) 


Shock Point Calculation 

In computing a shock point , the point k is first located geometrically 
on the shock wave, as in figure 4. In a region in which the upstream condi- 
tions are free stream, this point is found simply by adding a constant to the 
x coordinate of the previous shock point. In a region of variable upstream 
properties the shock point is located geometrically at the intersection of the 
shock wave and the opposite family characteristic line in the upstream region. 
The upstream properties are then found simply by linear interpolation. It 
should be noted that characteristics in the upstream region are dropped when 
they are no longer required in the solution (such as line AB in fig. 4). Some 
computing time is thereby saved. 

In order to find the properties of the shock point, k, another point k 
is computed as the intersection of the shock wave and the characteristic line 
from the point k - 1 (see fig. 5 )- This intersection point is determined by 


substituting the shock -wave angle 6 for p and the upstream region stream 
angle ' S u for 6 in equations (2) and ( 3 ). For very weak shock waves , the 
angle p ± S is almost the same as the shock -wave angle and an intersection 
point k can occur upstream of the previous shock point or extremely far down- 
stream. In this case the shock point properties at k are set equal to the 
upstream properties of the point. A regular shock-wave point will be computed 
if the intersection point occurs downstream of the previous shock wave. 

The next step in_the procedure for computing the shock point normally is 
to find a new point k 1 . A line is constructed parallel to the first -family 
Mach line through the points k and k - 1. The line passes through the point 
k. The point of intersection of this line with the second-family Mach line 
through the points k' + 1 and _k - 1 is the new k' point for the shock cal- 
culation and is designated by k' in figure 5 • The properties at the point 
k' are formed by linear interpolation between the points k’ +1 and k - 1 . 

The computation for the shock-wave point then proceeds by an iterative 
solution. The equations involved are equations (l4) through (l 6 ) which were 
obtained from reference ( 4 ), and equation ( 17 ). 


sin 6 0 + b sin 4 0 + c sin 2 9 + d = 0 
- 7 sin 2 & 




b = 


M u 2 + 2 


Mi 


•u 


c = 


d = 


2M U 2 + 1 

Mu 4 

cos 2 & 


+ 


(7 + I) 7 - 1 

h M u 2 


sin 2 & 


M- 


■u 


U4J 


Sk = S u + 


7 - 1 \ 


In 


2/M u sin Q - (7 
7 + 1 


- D1 

— 1 

j 


-7 Zn 


(7 + 1)M U 2 sin 2 9 
L (7 - l)Mu 2 sin 2 0+2 


T.T TT 


4(M U 2 sin 2 9 - 1 )( 7 M U 2 sin 2 0+1) 


'u 


(7 + 1) 2 M U 4 sin 2 0 


S k = & k ' ± A k ,(W k - W k .) + C k , ± D k ,(S k - S k .) 


(15) 


U6) 


(17) 
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As a first approximation to 5k, equation (17) is computed assuming the terns 
involving W and S are zero. A value of sin 2 6 is then obtained from Equa- 
tions (Ik) by solving the cubic by a trigonometric method. This trigonometric 
method is described in reference 5* Of the three roots obtained, the middle 
root is the desired weak shock solution. Equations (15) an d (l6) are then 
used to compute new values for Sk and Wk« With these values a new value of 
6k may be obtained from equation (17). If this value of 6k agrees well with 
the previous 6k confuted the iteration is terminated. If it does not, the 
new value of &k is averaged with the previous &k and the iteration con- 
tinues by recomputing sin 2 0, S k ,and Wp as before. 


Body Point Calculation 


The first step in computing a body point is to locate it geometrically as 
the intersection of a characteristic line and the body. The body may be either 
the centerbody or the cowl. If the body is supplied in tabular form a search 
is performed in the table and the location is found by the simultaneous solu- 
tion of two linear equations in x and r and the stream angle &k is obtained 
from the table. When the body is supplied in the form of a function, the loca- 
tion of the point is found by combining the equation for the body, equa- 
tion (l8), and the characteristic line equation, equation (19), to form equa- 
tion (20): 


r k = f ( x k) ( l8 ) 

r k = r k’+i+( x k - x k , +1 )tan(p + &) k , +1 (19) 

g(x k ) = f(x k ) - r k , +1 ± (x k - x k , +i )tan(q + &) k , +1 = 0 (20) 


Equation (19) may then be solved for by means of the Newton -Raphs on tech- 
nique (ref. 6). The stream angle is set equal to the arc -tangent of the 

slope of the body at the intersection point. 

The remaining properties are easily found since the entropy S^ remains 
constant on the surface between shock wave impingements. Thus, given 5^ 
and may be found from the compatibility equation, equations (5)- Equa- 

tions (8) through (13) furnish the remainder of the calculation. 


Starting the Solution 

The calculation may be initiated in two ways. The first method is to 
approximate a conical flow at the nose. The second method is to use some other 
procedure to calculate several points along a vertical input ray. 


8 



In the first case, the stream angle and the Mach number on the body and 
at the shock must be supplied to the program. From these quantities the 
remaining properties may be found. As few as two points may be used to start 
the flow field. However, three points are better when the Mach number is low 
(bejow Mach 2) or the centerbody angle is small (e.g., 9c = 5 °) • The third 
point is the average of the first two. 

In the latter case, any method may be used to compute a vertical ray 
consisting of a number of points at any station on the centerbody ahead of the 
cowl lip. The quantities, entropy, S^, velocity, W^, and stream angle, & k „ at 
each point are the required input parameters. 


Starting a New Region 


A new region is started after a shock wave has impinged on a body or after 
the flow field in the initial region has intersected the cowl lip. A two- 
dimensional flow is assumed in the immediate area of the intersection for the 
purpose of computing the starting ray. The initial ray has two points, a 
shock point and a body point both with the same x coordinate. The x coor- 
dinate is determined by adding a constant to the x coordinate of the inter- 
section point. This constant is determined by taking a percentage of the 
constant used in starting the previous region. The r coordinate at the body 
is determined from the body equation or table, as is the slope of the body. 

The slope of the body determines the stream angle at both points . The shock 
wave angle is computed by equations (lh). The entropy and the velocity may 
then be computed from equations (15) and (l 6). The r coordinate of the shock 


point 


H C or 


)Utcd using "til 6 SiiOu A.— wB,"V 6 Ullgre • 


A means of controlling accuracy in the program is by modifying the dis- 
tance from the cowl lip to the initial ray in the second region. This distance 
is computed as a percentage of the bow shock wave spacing. The percentage 
used is an input quantity. However , if none is provided, a value of 0.50 is 
used. The accuracy of the solution has been improved by use of a spacing com- 
patible with the local mesh size. 


Control of the Mesh Size 

Control of the mesh spacing has been found to be desirable in order to 
maintain a good distribution of output at the throat of the inlet. The char- 
acteristic mesh size in the initial region is controlled by the spacing on the 
input line and by spacing the shock points evenly along the bow shock wave. 

The latter technique is employed to limit spreading of the mesh. The spacing 
along the bow shock wave is controlled by an input quantity. The horizontal 
distance along the shock wave is computed as the product of this input quan- 
tity and the distance from the nose of the centerbody to the initial point on 
the input ray. 

The characteristic mesh in subsequent regions is controlled by the way in 
which the shock wave points are computed. The points are located geometrically 
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as the point of intersection of a characteristic line in the upstream region 
and the shock wave. The mesh size is thus controlled by the mesh size of the 
previous region. 

A means has been provided for discarding intermediate points in a verti- 
cal ray when it is desired to expand the mesh rapidly. This expanding pro- 
vides for increasing the accuracy in the nose region without substantially 
decreasing the speed of the solution. The discarding of points is under the 
control of the user by input quantities. 

Refining the Solution 

When the distance between points along a characteristic line exceeds a 
value specified by the user, a means has been provided for improving the 
accuracy of the solution as in reference 3 without increasing the number of 
points in the mesh. The values of x, r, S, W, and & at the points k and k 1 
are averaged together and the k 1 point is replaced by the averaged values. 

The same is done with the k 1 + 1 point. The k point is then recomputed. If 
the new value of Wp agrees closely with the previous value of the itera- 

tion is complete. If not, the points are again averaged and the point k 
recomputed. This is done until the value of W k does not change appreciably. 


Coalescence 

The coalescence of characteristic lines to form a shock wave is a serious 
problem, for it indicates that the inlet design is not a very good one. The 
user must then redesign the body contours and recalculate the flow field until 
the coalescence has been avoided. 

When two characteristic lines coalesce, the user is notified and a test is 
made in the program to determine whether or not the two lines have actually 
crossed. If they have, the downstream characteristic line is dropped, a mes- 
sage is printed, and the flow-field calculation continues. The user can then 
determine if the shock wave is in fact building up or if the coalescence 
occurred because of inaccuracies arising from some of the approximations used 
in the program. 

When two characteristic lines become arbitrarily close to each other, the 
accuracy of the solution tends to decrease. There is a significant loss in 
the number of digits of accuracy. To avoid this problem, it is assumed, in 
this case that coalescence has actually occurred and the downstream character- 
istic line is dropped. 


CONCLUDING REMARKS 

A Fortran IV computer program to calculate the flow fields in three- 
dimensional axisymmetric or two-dimensional inlets has been written. 

In writing the program, every effort has been made to keep the computa- 
tion as rapid as possible. The calculation of the upstream and downstream 
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regions has been carried on concurrently to avoid the computation of extra- 
neous flow-field points and the necessity of saving all the points in the 
upstream region. In this way the use of external storage devices has been 
avoided . Extensive search techniques and elaborate curve -fitting schemes have 
been avoided as well as high-order interpolation formulas. As a further aid, 
several cases may be stacked to be run at the same time. 


The program is flexible since the flow field may be initiated in two ways 
and the surface input may be tabular or analytical. When tabular input is 
used the compression surfaces can be defined by both coordinates and local sur- 
face angles. The input of surface angles is redundant j however, coordinates 
alone cannot usually be calculated for input accurately enough or smoothly 
enough to give a uniform mesh of characteristics. It is suggested that the 
surface angles be plotted and faired before they are used as input. 

The rapidity with which a case may be computed (from 0-5 to 2.0 minutes 
per case) has been an advantage in designing inlets, since a trial and error 
design procedure involving many cases is often necessary. 

The suitability and accuracy of the program for the design of supersonic 
inlets is illustrated by the inlet design shown in figure 6. The input and 
output for this case are presented in appendix A. This particular case was 
computed in less than one minute on an IBM 7094 computer. 


Ames Research Center 

National Aspon&u/fcicG 
Moffett Field, 


and Space Administration 
Calif., Mar. 25, 1965 
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APPENDIX A 


PROGRAM USAGE 

In this appendix, the program usage is detailed. An attempt has been 
made to keep the usage of the program simple. The program itself is listed 
in figure 7. A sample case is listed as figure 8. The input cards used to 
obtain the data in figure 8 are listed in figure 9* Figure 10 contains a 
listing of a sample plotting program and figure 11 contains listings of sample 
body contour programs. 


DECK MAKE-UP 


Deck 

label 

Subroutine 

name 

Description 

EF3131 

(main) 

Control program. 

EF3132 

BODY 

Computes a body point. 

EF3133 

FLOW 

Computes a flow-field point. 

EF3134 

JUGGLE 

Refines solution when mesh size is large, by averag 
ing upstream points. 

EF3135 

PUNT 

Output subroutine. 

EF3136 

CBODY 

Dummy subroutine. If the centerbody ordinates are 


given in analytic form, this program should be 
replaced by a subprogram that computes the function 
and its first derivative. The function should be 
in the form, r = f(x), where r is the distance 
along a line perpendicular to the body axis and x 
is the distance along the body axis. The CALL 
statement is "CALL CBODY (l,X,R,DR)" where X = x, 

R = r, and DR = dr/dx. At the beginning of the 
control program an entry is made to this subprogram 
with the control word I = 1 to allow for the ini- 
tialization of the subprogram. For all other 
entries 1=2. On the initial entry X should be 
set equal to the last value of X on the center - 
body. 

EF3137 ABODY Dummy subroutine. If the cowl is given by an ana- 

lytic expression, this program should be replaced 
by a subprogram that computes the function as 
indicated under CBODY- In this case the CALL state- 
ment is "CALL ABODY (l,X,R,DR)." On the initial 
entry, when 1=1, the cowl lip ordinates should 
be returned in X and R. 
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Deck 
label * 

Subroutine 

name 

Description 

EF3138 

k 

SHOCK 

Computes a shock point. 

EF3139 

CONIC 

Sets, up a conical input ray. 

EF3140 

ERROR 

When an error condition is encountered, this sub- 
program prints an error code which may be used to 
determine the type of failure. The next case is 
then read in under most circumstances . If the error 
is irrevocable, EXIT is usually called. 

EF3141 

ENDFIL 

This subroutine is entered when the last stacked 
case has been completed. This condition is sig- 
nalled by the input card with the word "DONE" in 
columns 1-4. 

EF3142 

FLINT 

Computes with the point of intersection of the flow 
field with the cowl lip. 

EF3143 

SPLOT 

Dummy subroutine. To be used if plotting is not 
wanted. (See EF3150.) 

EF3144 

ACRAY 

Computes a two-dimensional input ray. 

EF3145 

CUBIC 

Finds roots of a cubic and selects proper root for 
the shock angle equation. 

EF3146 

UPSC 

Computes the upstream conditions for a shock point. 

EF3147 

BSINT 

Computes the intersection of a shock wave with 
either the centerbody or cowl. 

EF3148 

CBODY 

Sample subroutine for computing the centerbody ana- 
lytically. 

EF3149 

ABODY 

Sample subroutine for computing the cowl analyti- 
cally. 

EF3150 

SPLOT 

Constructs arrays for plotting. The user must 
supply a program that writes the plot tape. 

EF3152 

PAGE, TITLE 

Controls page numbering and headings . 
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TAPES 


Logical tape no. Usage 


5 INPUT 

6 OUTPUT 

7 PLOTTING OUTPUT (Optional) 


INPUT CAPDS 


All input is in floating-point form, 7 values per card, 10 columns per 
value, except for the first card. (See fig. 9 Tor a sample case.) 

Card no. Columns Contents 


1 

2 


1-72 Case identification in Hollerith. 

1-10 N = number of input points if initial ray is to be read 
in. N = 0, if initial input ray is to be computed as 
conical input ray. (N < 50) 

11-20 NOB = number of centerbody ordinate points if body is in 
tabular form. NOB = 0, if body is analytic. (NOB < 50) 

21-30 N0A = number of cowl ordinate points if cowl is in tab- 
ular form. N0A = 0, if cowl is analytic. (N0A < 50) 

31 -L0 MAXNP = maximum number of points allowed in a ray. If 

no value is specified, it is set to 50. (MAXNP < 50) 

1,1-50 NOBIS = number of points to be discarded between the 
points kept. When the number of points in a vertical 
ray exceeds the value assigned to MAXNP, and N0DIS = 1, 
every other point is discarded. 

51-60 IDU'l = 2 or 3 determines whevher the solution is to be 
two-dimensional or three-dimensional axisymmetric . If 
no value is assigned, a three-dimensional solution is 

cl 3 3 LU.'iCCL * 


1-10 EMIN, M„. - free- 3 trearri Mach number. 

J \J*j 

11-20 GAMMA, 7 = free -stream ratio of specific heats. 

21-j0 SING, S - free- stream entropy, dimensionless . 


Ik 



Card no. 


5 a 


5b 

6 

7 


8a 


8b 


y 

10 


Columns 

Contents 

31 -bo 

THETA, = shock wave angle (in degrees' 

1-10 

TEST, convergence test for iterations (l0 _fc gives good 
results). 

11-20 

CRMAX, maximum distance between mesh points. If this is 
exceeded, a refinement is attempted. 

21-30 

COALT, coalescence is said to occur if the increments in 
x and r are both less than this test quantity. 

31 -b 0 

SPACE, controls spacing on the shock wave (in the x 
direction) , as a percentage of the distance from the 
nose to the initial input ray. 

bl -50 

SPACC, controls distance from the cowl lip to the initial 
ray in the second region as a percentage of the bow shock 
spacing. If no value is given SPACC is assumed to be 
0 . 50 . 


(Optional - include if body is defined in tabular fashion.) 
x ordinates - NOB values on the centerbody. 

(Optional - may be included or not if centerbody is defined by an 
analytic function. ) 

User's format - Constants used in centerbody program. 

(Optional - include if body is defined in tabular fashion.) 
r ordinates - NOB values on the centerbody. 

(Optional - include if body is defined in tabular fashion. ) 

S angle (in degrees) that the line segment joining adjacent 
table values makes with the body center line. NB = NOB-1 values 
should be specified. 

(Optional - include if body is defined in tabular fashion.) 
x ordinates - NOA values on the cowl. 

(Optional - may be included or not if cowl is defined by an ana- 
lytic function.) 

User's format - Constants used in cowl program. 

(Optional - include if body is defined in tabular fashion.) 
r ordinates - NOA values on the cowl. 

(Optional - include if body is defined in tabular fashion.) 

S angle (in degrees) that the line segment joining adjacent 
table values makes with the body center line. NA = N0A-1 values 
should be specified. 
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Card no. 


Columns 


Contents 


11 


11 


12 


(Optional - include if program is to compute the initial conical 
input ray, that is, N = 0.) , 

1-10 XB, x Q , initial x value on the Body at which the mesh 
is to start. 

11-20 RB, r Q , initial r value on the body. 

— - — — - - ... i. . _ - , — v. -I- Vi i— v V T.rPl ~\T . 

21-30 EMbtt, W 5 H ocv-’ lvia ' uli lluiliucJ - W11 

31-40 EMBODY, Mach nuiriber on the ’ bod Y- 

41-50 DELSH, & SH0CK degrees). 

51-60 DELBD , & B0DY (in degrees). 

61-70 N0PIN, number of points in initial conical input ray. 

(Optional - include N of these cards if the initial ray is to 
be read in, instead of being computed by the conical input sub- 
routine . ) 


1-10 


local 

X 

11-20 


local 

r 

21-30 


local 

& (in degrees) 

31-40 


local 

W 

41-50 

s j ,*■’ 

local 

s 

1-4 

The word "DONE." Any number of cases may be 


before this card. This card is needed on the last case 
in order to complete the plotting arrays and terminate 
properly. 


INPUT FOR PLOTTING 


These cards may be in any format the user chooses. The sample program 
included uses a 7 FIO .6 format and is set to read 5 cards placed after card 
no. 1 and before card no. 2. These cards contain the origin of the plot, the 
scale factors, the maximum value for each variable, and the minimum value of 
the ordinate. The Mach number and pressure distribution along both bodies are 
plotted versus x. The mesh is also plotted. 


16 



ERRORS 


When an error occurs during the computation of a case, an error message 
is usually printed and the program goes on to the next case. 


The 

following is 

a list 

of error messages that may appear: 




Error Conditions 

Deck 

name 

Subroutine 

name 

Error 

code 

Probable cause of error 

EF3131 

(main) 

1 

The number of points in a ray exceeds MN (an 
input quantity < 50) or 50 if MAXNP is < 0 . 

EF3132 

BODY 

2 

The iteration scheme for finding the intersec- 
tion of a Mach line and the body is not working. 
The subroutine for computing the function and 
its derivative may be incorrect. 

EF3132 

BODY 

3 

The local stream angle, &, is greater than jt/ 2 . 

EF3132 

BODY 

1 

M* 

The iteration for a body point is not converging. 

EF3133 

FLOW 

5 

The iteration for the local stream angle, 5 , is 
not converging . 

EF3135 

PUNT 

6 

There is an error in the computation of the 
local Mach number. 

EF3135 

PUNT 

7 

The recovery, P t /P toQ , is greater than 1 . 0 . 

EF3138 

SHOCK 

8 

There is an error in computing local entropy, 

S • 

EF3138 

SHOCK 

9 

There is an error in computing local velocity, 

W. 

EF3138 

SHOCK 

10 

There is an error in computing the shock wave 
angle, 9 . 

EF3139 

CONIC 

11 

There is an error in computing the entropy, S. 

1.1 1 .1 ^ -1 \ 

ShU $±*+<1 

FLINT 

12 

The shock wave has fallen inside the lip. 

EF 3144 

ACRAY 

13 

There is an error in computing the local entropy, 

S. 
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Deck 

name 

EF3144 

EF3144 

EF3146 

EF3146 

EF3147 


Subroutine Error 
name code 


Probable cause of error 


ACRAY 14 

ACRAY 15 

UPSC 1 6 

upsc 17 

BSINT 18 


There is an error in computing local velocity, 

W. 4 

There is an error in the shock wave angle, 0, 
or the Mach number, M, in the computation of 
a new region. The upstream Mach number is less 
than the local Mach number or the local Mach 
number is subsonic. 

There are too few points in the upstream region. 

There is an error in the computation of the 
upstream Mach number. 

There is an error in finding the point of 
intersection of the body and the shock wave. 
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APPENDIX B 


FLOW DIAGRAMS 


A. complete set of flow diagrams is included in figure 12 as a means of 
assisting the program user . These diagrams were drawn to conform to standard 
flowcharting techniques (ref. 7). Each diagram has been identified by the 
corresponding subroutine name and associated deck name. 
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Figure 5*- Location of field point used in the shock point computation. 



Figure 6.- Typical inlet design. 
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Figure 7*- Program listing. 
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Figure 7 • * Continued 
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Figure 7 • - Continued 
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Figure 7 . - Continued 
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Figure 7*- Continued. 
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Figure 7-- Continued. 
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Figure 7 . - Continued 
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Figure 7*- Continued. 
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Figure 7.- Continued- 
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Figure 7 • - Continued 
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Figure 7.- Continued. 
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Figure 7.- Continued- 
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Figure 7.- Continued 
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THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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Figure 8.- Sample output • 
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Figure 8.- Continued 
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Figure 8.- Continued 
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Figure 8 . - Continued 
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Figure 8. - Continued 


THREE-DIHENS ICJNAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 

CASE 14A 12.5 DECREE CONE 
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Figure 8.- Continued. 



THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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'igure 8.- Continued. 
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Figure 8 . - Continued . 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE. 03/12/65 
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Figure 8.- Continued. 
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Figure 8.- Continued. 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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Figure 8.- Continued. 


three-dimensional axisymmetric method of characteristics, internal case 
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Figure 8.- Continued. 
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Figure 8.- Continued. 


THREE-DIMENSIONAL AXISVMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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Figure c . - Continued 
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Figure 8. - Continued. 


three-dimensional axisymmetric method of characteristics, internal case 
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Figure c . - Continued 



THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 

CASE 1 AA 12.5 DECREE CONE 
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Figure 8.- Continued. 


•J- -HoaJi^^Mifun^iAaoovoj-^HO 

z^cOrHMf\jo-iintoo^f\jmfriir, o^i^i 
•-■ON^ i^inin>o^^Nr-^SNNso(c 
hJ’^pawfOffiro^fnmmrfwnmmNcoN 

Q.O' 5 ' 0 » 0 '(M>(MM^a l (M>oo 5 tC offl ffl 

N 0 ' 0 (>i 7 ' 0 'a' 0 ' 0 'a'( 3 '(j'(j' 0 ot>{> 5 '^ 


S££££ r ^'°'£ 0 , ^ (NJinornirv( >CF'®®»'-ir' 

^^©^^ininmsONO^NssNNoooovO 

O'CT'O'^ooooooa'oooCT'a'a'oa'oo 


0.00000000000 000 00 0 0 ooooooooooooooooooooo 


• OfflNoO^rftaoONON'NJOONftiflSaJ 

amo^»Njo>p^o<Mh-oocn^tnir\r^rMh-m 

^OO»-i»^^rMr 0 m^tn^s 0 r-a 000 crN<r»fO 

U>0'0'0'0«0<fl'0'0'0^>0^0v0'0^>j-J- 

5**** 1 ••••••••••••• 

X:rgr\jfMf\i(\j<Njrgrgrsifv»ro«Mru(MCMcgr>jcM 


ooirvLf\^-i^'mmcor«.r^(\j'^fMrgrriQOfsjirvsO'< 0 ^ 

O'H^MMoomHinO'^coMtisHO^s^ 

!T ^ 5 ? N '°" ,fWf,, ' 0 ' 0 rnHco ' 0 't^' 0 'a^ 

'T^ -T ' 0 'OvOO^vO>OvO^CvOOO>OJ-^»)-J- 

••••••••••• 

rsjfMfNirg<\jf\jf\j(\jrsjrMrgfgrsj<\4(Mc\j<M<MCMfM(M 


O 

o 


r^S^^ LA *> ,0:7 ' rWs0 '^ Ln ' orNJ ' 00rs -^ 1 ^ 

ULJO(\j<>rsiNnNO'j*No<OrOiAfOfflN<t 

QO-tO-Offl'Oins^OOMt^CDtT'N^vO 

-OHiONNON^-NN'O’Off'-i-^Om® 

<OrHOm(\jOrrirri v o^HQor«-h- 3 'r\jor^fM 

h;lABHOHv 0 (\ J 00 s J. P HS^HC 0 O(>HN 

• • • 


^*^'^I^?r' 5 rsJr ^ c0Oc0 ' , '' 0 '-'' n ' 0 ^^o 

N^ins^socooj^oomoj^mO'OOs® 

Lnr04>rvjf\jh-rrifMa3C0c\jr»-sd-{j%-o(M>00'CD(Nj 

'J'^'iAOmOrNjr-sDcomOh-f^CD^HONDOjD 

^a)HlA 3 '^ 7 >^ 0 ^rl^^ rnO( j ) OHlA 0 


dddd 2 dddddddl^hV*I,o 2 dd- 


z 

* 


W®{MM^ 4 ’'OvOoOO'HMfloO'OOi/\(M 
ILNinoJ-NH^WxflflONOJlrjrgNON 
ZO'-(sOCT'fNj^Ofsi^irkO>o^»>i'C 7 'fsj ®^4 
— 'W 3 'Oooin(j'fOvo©^OH(\j(\j(vjOMn 
aHOoooMn^MOooNirtfnHa'^om 
'v.coaccor-h-r-r»-r--r > ->OsOv£)£>vOunrornm 


vi-vi-crrriOh-ir>fNjfN(<Njr^r«.ofNj— i^^r^-tsOr -4 
®CDfl<)CO^in« 1 ,i-SOi\jO«OTl! 7 '^HOJ-N 
nN 3 ' 0 'Mma'NH^fNj f r l o>rO^OOJ' 0 ^jaD 
^nO'Hin^ogcor<ij)a'H(\^invOv 0 ^oo'(n 
C7'a'inMO!>GOvOlCfr|HOQOsOvt(\JOCOrNjrJCD 
mmmaocoKh.rwf«-h-r**h->o«Osa' 0 >orn^^-m 




|s||ps|£gs5l5g|sssI”sssssSSl35Hg"lSc223 

«sss2SsSElsIlsHiss?ipss;si2'^ 


o o o o 


. . ^ trCrOCr ' a ' Ln ‘ A ' 0 <> ' 0 r '^^®®®^ 0 ^^' 0 CT' 0 ' 

Dooooooododddddd^dddddddddddddddddddd 


Q 

a 

X 


22 I\i 22 J 21 f\!Si'S^ 2 ’ 5 « ,a ' ,i *’ ro<vl ' d ’ , ^ a ' 00 <Ni 00 m ‘ n ' org '- ,orr ' ( ^'t< 7 '(y>oo'CT' 

S*Ssgss5SS*llll|liIsll|isS5Sssfis5S2KSSS 

K 'tt’V r .t^* 23 JS 5 SS?S 5555 ?i 5 S 55 J JSSSSSSKSSSS 


)r \)i‘^'tir\\ONCO' 7 'Of-<f\Jfr\^ifi^ 




X 

< 


<' 0 - 0 > 0 ' 0 ' 0 ' 0 > 0 ' 0 >OnO' 0 > 0 ' 0 ' 0>0 


r^r-^o<)' 00 ' 0 s 0 o>o^' 0 > 0 ' 0 v 0 v£>h-cocooc 

'O'O'O'OnOsOsO^nOnOOnOO-O 


l N >\l ^J fNI r\j r\| r\j ^ ^ , 


Q 

t 

Ui 

UJ 

QC 

X 


81 


Figure 8.- Continued 



three-dimensional axisymmetric method of characteristics, internal case 

CASE 14A 12.5 DEGREE CONE 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8 . - Continued . 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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Figure c . - Continued 



THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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Figure 8.- Continued. 
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Figure 8.- Continued. 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS, INTERNAL CASE. 
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Figure 8.- Continued. 
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Figure 8.- Continued 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued 
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Figure 8.- Continued 


THREE-DIMENSIONAL. AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE. 03/12/65 
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Figure 8 . - Continued . 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE. 
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Figure 8.- Continued 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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Figure 8.- Continued 


THREE-D1MENSI ONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE. 
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Figure 8.- Continued. 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS, INTERNAL CASE. 
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Figure 8.- Continued 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE. 03/12/65 


[ 


2 r 

UJ 

o 

on 


< 

o 

o 


u.rgNa'iM^oooNOirivOs^^tnOtc^ojoinsOrnm^ONao 

>ia5inOO'(MNNN^rtNmM»HCOinfOin^(nN(M7'CO>J , f-*0 
kH4 , 00'Ol/MnOOOKO'H(f|>ON'0<t^ONS^rJin'OlAOJHCf l 
QLNp*-'0Nr>»Ni* , -‘0'0'0r , “f , ''0r , -^NNr'>0\0ON >o s n n 
s5'!T , (J'C>(J'0!J'0'13'0'!7'C>0'0'C'!>0'Off't7'0!T' ^ O' O O' O O' 


>0 00 N O' m ^ (fl i-l H vO o 

oinrgHNooff'MnoH 
co-Oif— lO'-grop-co-im 
S'OKiNO'omcoooN 
N N S N ifl >0 O ■O S >0 N 
^O'^O'O'^O'O'CP^ ^ 


CLOOOOOOOOOOOOOOOOOOOOOOOOOOOd 


o o o o o 


o o o o o 


o 


• 

if lO 

if o 

m 

O' in n 

sO P- 


p-l P~ 

ao 

ao >f 

.-4 

O' CO 

o 

o if 

-h e- 

r- 

rO O' «o 

if if 

O' 

O' P- 

if 

— i if 

in 

p-4 p- 

z 

if r» 

in cc 

n 

in in co 

r -4 O 

o 

•f P-4 

o 

o oo 

eg 

cr ro 


iO if 

in «f 

m 

f Htn 

O' O' 

r- 

eg csj 

iO 

in if 

o 

co eg 

X 

oo f- 

*h nj 

Csl 

— * -O —i 

if if 

oo 

oo eg 

CNJ 

rg O' 

o 

eg eg 

a 

in in 

ro ro 

ro 

m rfi rn 

in in 

in 

m ro 

ro 

ro eg 

ro 

ro ro 

X 

■— i p— i 

H r-i 

-4 

H P —4 pi^ 

H H 

p — 4 

H H 

p-i 

1-4 —4 


■—4 p— * 


iO 

r— .-I 

—4 

n ao 

iO 

r0 

eg 

P- 

O 

ro 

O' 

O 

O' 

in 

eg 

— 1 

o 

eg 

if 

■f O' 

if 

00 o 

>0 

in 

in 

p- 

iO 

ro 

p—4 

Ph 

iO 

O' 

pH 

-—4 

n 

ro 

p-4 

o o 

p- 

•H if 

in 

CO 

if 

in 

pi- 

P- 

o 

ro 

o 

n 

pH 



ro 

O' 

O' f- 

H 

lO lO 

iO 

m 

in 

pi* 

pi- 

P- 

if 

O 

n 

eg 

n 

in 

o 

O' 

if 

if ao 

eg 

eg o 

if 

ro 

rO 

eg 

O' 

if 

p—4 

>f 

ro 

co 

if 

if 

O' 

O' 

n 

n n 

ro 

ro nj 

ro 

ro 

ro 

ro 

l\J 

ro 

ro 

ro 

ro 

in 

in 

in 

eg 

eg 

—4 

p— 4 p— 4 

pH 

•H p-4 


—4 

J, 

p-4 

-H 

—4 

JJ 

J, 



p—4 



^■i 


Jom!7 i in^-oiniM^iftcoffi(MON(n-(fn o >o 
UI'OCOON'O'tfOOtMfl^COOOtM^fl^NJ' 
OooONcoNfMin^cO'OHNomrjtn^in^ 

— ^■OiOJ-NB>0>J->}-Oinco^O'O^^HS^ 

< J- O co co e ifOtNjONcofloos-t^in^o 

hNH^OMOffl^(nO(MmNOHin>O^Mn 

O uj^^oooooonoocooO'-'hcocooonoonooo 

o o I I I I I | | | I I I I I I I I I 


in oo 

if oo o >o o n 

3 rg 

iO 

co O' 

in 

O ro Ph 

o 

ro 

n ro 

O' O in O' eg ro 

O O 

o 

O' 0 - 


O O eg g 3 

eg 

co o 

O' eg ro ro ro rg 

O pH 

O' 

iO eg 

iO 

O -< ro 

co 

eg 

O' n 

ro O O P— e- iO 

O if 

O' 

n ® 


O >0 n 

PH 

n 

p- ao 

iO O' —4 ro n -o 

O eg 

eg 

in iO 

if 

O >o in 

oo 

ro 

O rg 

co o if ->< pH eg 

O ro 

O 

eg o 

ro 

pH P-4 O 

eg 

eg 

O r -4 

ao oo rn o' pH co 

ao e- 

O' 

ao co 

00 

0 . 

0 . 

0 . 

7 , 

PH 


* I I I I I I I I I I I I I III 


M 

X 


UJ 

z 

3 

o 


UJ 

UJ 

O' 

CJ 3 

UJ 


Z 

CL 


oo i n 

O' o 
o if 

in e* 


ro m 


if OOO'Of 

O eg if 

eg 

p-4 oo eg id m o 

p-4 m 

eg 

i f eg 

oo ao 4- 

>f 

if 

O 

in 

>f 

O' 

iO o 


eg go 

o rg 

eg 

O' pH if 00 if 

eg pH O 

eg 

co p—4 co rg in if 

•f O' 

PH 

iO eg 

co o co 

O' 

p—4 

•o 

ro 

co 

eg 

if p-4 

co 

fH 

co 

O' ro 

O' 

O' ro — 4 O' ro 

ph co o 

if 

rg ro iO eg eg O 

eg iO 

CO 

O if 

~4 iO ao 

1 0 

eg 

p-4 

■f 

ro 

co 

nO if 

00 

PH 

if 

O -4 

o 

if O' ph in if 

co ao o 

in 

eg if O' i0 if P- 

m ao 

CO 

o in 

-4 iO O' 

<f 

rg 

co 

ro 

fH 

'O 

ro O 

ro 

ro 

co 

O p—4 


o ^ p-t in eg 

iO O ro 

in 

ro ro O O —4 in 

iO ro 

o 

ro m 

rO O' O' 

O' 

O' 

ro 

CO 

co 

co 

if o 

ro 

gj 

lO 

O' co 

Ph 

in in in g) co 

in o o 

in 

o if in in o ph 

ro if 

o 

o in 

if f O' 

o 

eg 

eg 

<f 

O' 

o 

iO -4 

eg 

in o 

O O' 

O' 

eg eg eg rg .-4 

eg O' O' 

CO 

ao eg eg eg ro <-* 

p—4 —4 p—4 p— | p—4 

12. 

12. 

o> 

O' oo 

rg rg eg 

*H «H pH 

eg 

r*-i 

eg 

eg 

r*^ 

eg 

eg 

«-! 

eg 

.—4 

eg eg 

sH — < 

eg 

co 

O' 

O' eg 

eg 


Q 


in 

. 

eg 


< 

>f 


UJ 

on 

< 

o 


—4 

—4 

rg 

P -4 

o 

pH 

O' 

o 

ro 

in 

in 

if 

eg 

O' 

if 

lO 

ao 

co 

o 

ro 

m 

pH 

eg 

o 

eg 

eg 

-o 

Q 0 

in 

O' 

O' 

iO 

o 


o 

fH 

in 

if 

ro 

O' 

if 

oo 

ro 

in 

if 


o 

in 

rn 

eg 

o 

eg 

m 

i 0 

pH 

O 

ro 

o 

in 

h* 

ro 

o 

pH 

PH 

ao 

ro 

if 

00 

o 

o 

pH 

CO 

if 

o 

pH 

PH 

if 

if 

-o 

if 

p -4 

ro 

>0 

ro 

pH 

o 

if 

4 f 

ao 

ro 

o 

eg 

eg 

O' 

|H 

»f 

g 3 

if 

if 

ro 

pH 

0 

eg 

iO 

O 

m 

co 

CO 

O' 

ro 

O' 

eg 

o 

o 

o 

ro 

O' 

oo 

iO 

O' 

rg 

O' 

O' 

eg 

4 f 

eg 

O' 

co 

in 

iO 

o 

co 

m 

if 

O' 

if 

eg 

O' 

ao 

PH 

■o 

ro 


eg 

*— 4 

ro 

ro 

o 

O' 

ao 

if 

in 

GO 

iO 

if 

o 

if 

ro 

pH 

o 

oo 

in 

if 

oo 

O 

pH 

ro 

rg 

eg 

H 

O' 

iO 

in 

eg 

00 

O 

pH 

#H 

pH 

o 

fH 

iO 

if 

eg 

rg 

pH 

ao 

sO 

if 

ro 

p- 

O' 



eg 

CO 

GO 

O' 

O' 

oo 

co 

oo 

PH 

CO 

oo 

00 

ao 

O' 

O' 

CO 

co 

ao 

co 

PH 

co 

00 

CO 

O' 

O' 

co 

co 

co 

co 

O' 

O' 

00 

co 

co 

CO 

P- 

Ph 

co 

00 

O' 


• 


. 



• 

• 

. 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


. 

• 

* 

• 

• 

• 


. 

• 

* 



• 

« 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 


o 

O' 

rO 

iO 

eg 

O' 

in 

o 

4* 

in 

^-4 

O' 

® 

43 

O 

O' 

if 

O' 

in 

if 

in 

H 

in 

O 

>f 

eg 

ro 

in 

ro 

in 

>o 

if 

'f 

ro 

in 

iO 

co 

in 

-f- 

in 

o 

ro 

if 

eg 

pH 

(NJ 

eg 

PH 

O 

o 

PH 

PH 

UD 

4 ^ 

o 

pH 

H 

ro 

ro 

iO 

o 

PH 

O' 

eg 

in 

O' 

o 

4 - 

ro 

iO 

in 

co 


>0 

in 

if 

>0 

Ph 

ro 

ro 

03 

ro 

co 

m 

p-4 

ro 

fH 

rO 

if 

O' 

O' 

® 

if 

if 

ro 

PH 

o 

® 

ro 

if 

O' 

O' 

if 

rg 

<£> 

eg 

O' 

H 

eg 

o 

ro 

PH 

ro 

ro 

eg 

O' 


O' 

O' 

O' 

H 

PH 

H 

eg 

O' 

eg 

O' 

® 

iO 

lO 

® 

w 

H 

o 

CO 

p-4 

® 

rg 


o 

o 

ro 

fH 

iO 

O 

e- 

o 

H 


if 

ro 

>o 

O' 

O' 

iO 

PH 

ej 

rg 

o 

in 

>o 

o 

in 

P“ 

o 

i0 

iO 

m 

sf 

in 

iO 

PH 

iO 

® 

pH 

P- 

PH 

IH 

i0 

if 

in 

o 

lO 

ao 

ro 

CO 

ro 

co 

PH 

O' 

co 

O' 

ro 

ro 

ro 

in 

iO 

fH 

P- 

O' 

o 

o 

iO 

o 

g> 

NO 

lO 

>o 

o 

iO 

lO 

lO 

iO 

iO 

>0 

<3 

iO 

i0 

iO 

lO 

iO 

iO 

NO 

iO 

iO 

10 

<o 

iO 

lO 

iO 

iO 

<3 

ro 

ro 

iO 

iO 

lO 

10 

lO 

>o 

• 

• 


• 


• 

. 



• 


• 

. 

• 

44 


• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

. 

. 

• 

• 

• 


• 


• 

. 

• 

* 

• 


ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

rO 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 


►~iftng3.Hegroifin.Hrgro.f«o-Hogroifinro.-4rgror-rHpgmiftn>o.HegroiftniO«HrgrofHrH 

O 

Q. 


V O' O' 

<mm 


4 ‘ififif'fOOOO.f 

pHpHpHp— 4»— lifpfij-sji.— 4 


minininininoooiris0s0>0>0'J3*o 

pH H <— I p—4 •— I pH >f if if pH »— I pH i— I •—! H t— 4 


pH f- 00 

if pH pH 


f f f f m 

UJ 

ae. 




ro m ro 


>f if if 


if if if if if 


«f if if if m 


ra if if 


105 
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Figure 8.- Continued 
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Figure 8.- Continued 
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Figure 8.- Continued. 
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Figure 8.- Continued 
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Figure 8.- Continued 


THREE-DIMENSIONAL AXISYMMETRIC METHOD OF CHARACTERISTICS. INTERNAL CASE 
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Figure 10.- Sample plot program. 
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Figure 10.- Continued. 
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Figure 11.- Sample body program. 
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Process information or compute 




Connector 


Terminator 


(a) Symbols used in flow charts. 
Figure 12 . - Flow charts . 
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BODY 



(c) Deck EF3132 subroutine BODY. 
Figure 12.- Continued. 
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PUNT 




(f) Deck EF3135 subroutine PUNT. 
Figure 12.- Continued. 
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SHOCK 


Save 

previous 

9 


wave /body 
intersection. 


Replace shock 
point by point 
of intersection 



Interpolate to 
to find j'k' point 
for the shock point 


x Two \ 
dimensional 
\case/ 

iNo 


Compute 
T, by eq (4) 


Compute point 

x j k » r jk as 

intersection of shock 
wave and first 
family characteristic 


Compute 

upstream 

conditions 


Compute 
p, by eq (4) 
C|by eq (la) 


Compute 
•jr by eq (la) 
Dj by eq (la) 


Compute 
S| by eq (17) 


X^hock wavex 
and characteristic 


Set local 

properties equal 
to upstream 
properties 


Compute 
sin 2 0 by eq (14) 


Compute 
S j k by eq (15) 

W jk by eq (16) 


Compute 
St k by eq (17) 


Ax, on shock 
vwave <0 s 



Set spacing, Ax, 
on shock wave 


Compute shock 
point as previous 

8 I 8 jk 

rr\ 1 - t — ^ ^ 

at 50% of 
initial x 


shock point 
plus Ax 

^ i> 


(g) Deck EF3138 subroutine SHOCK. 


Figure 

12.- Continued. 


has conveged to 
























CONIC 



ERROR 



Print the 
error message 





(i) Deck EF3140 Subroutine ERROR. 


ENDFIL 



(j) Deck EF3141 subroutine ENDFIL. 


(h) Deck EF3139 subroutine CONIC. 

Figure 12 . - Continued 
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ACRAY 


Second 

region 
































(p) Deck EF3150 subroutine SPLOT. 
Figure 12.- Concluded. 


NASA-Langley, 1965 


A -207 k 


131 





